Stationary localized states due to nonlinear impurities described by the modified 

discrete nonlinear Schrodinger equation 
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The modified discrete nonlinear Schrodinger equation is used to study the formation of stationary 
localized states in a one-dimensional lattice with a single impurity and an asymmetric dimer impu- 
rity. A periodically modulated and a perfectly nonlinear chain is also considered. Phase diagrams 
of localized states for all systems are presented. From the mean square displacement calculation, it 
is found that all states are not localized even though the system comprises random nonlinear site 
energies. Stability of the states is discussed. 
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I. INTRODUCTION 

Localized states appear due to the presence of impurity 
or disorder (which breaks the translational symmetry) in 
the system [0. While there have been studies on the 
formation of localized states due to linear impurities in 
various systems Q , only a few look into the formation of 
localized states due to nonlinear impurities. The discrete 
nonlinear Schrodinger equation [2-16] employed to study 
the formation of stationary localized (SL) states is given 
by 



,dC n 

! — ; — 

dt 



e n C n + V(C n+1 + C„_x) - X n\C n \ a C n , (1) 



where a is a constant, C n is the probability amplitude 
of the particle (exciton) to be at the site n, e„ and \ n 
are, respectively, the static site-energy and the nonlinear 
strength at site n, and V is the nearest neighbor hop- 
ping element. The nonlinear term ICn^Cn arises due to 
the interaction of the exciton with the lattice vibration 
||l5| , fl6|| . If the oscillators are purely harmonic in nature, 
(7 = 2. The Eq. (1), has been used to study the forma- 
tion of SL states in a one-dimensional chain as well as 
in the Cayley tree with single and dimer nonlinear impu- 
rities JlCHiai. For the case of a perfectly nonlinear 
chain where Xn = X f°r ai l n it has been shown that SL 
states are possible even though the translational symme- 
try of the system is preserved [j^Jl^-|l4| . These results 
were remarkably different from those of the correspond- 
ing systems with linear impurities. The Eq. (1) has been 
derived with the assumption that the lattice oscillators 
in the system are local and their motions are much faster 
than the exciton motion. A natural question to be raised 
is what happens to the formation of SL states when the 
oscillators in the lattice are harmonic in nature and cou- 
pled with their nearest neighbors. In such system, in the 
presence of the interaction of the exciton with the lattice 
vibration, the time evolution equation governing the dy- 
namics of the exciton has been derived by Kopidakis et 
al. as 

dC 

l ~df = ^c n + v(c n+1 + c n ^) 



Xn(\C n+1 \ 2 + \C n - 1 \ 2 + 2\C n \ 2 )C ni 



(2) 



where C n , e n , V and %„ carries the same meanings as 
in Eq. (1). The Eq. (2) is called the 'modified' discrete 
nonlinear Schrodinger equation. The Hamiltonian which 
generates Eq. (2) is given by 



H — (C*C n+ i + C n C* +1 ) 



n— — oo 



E f L (|c n+ ii 2 + |c„-ii 2 + 2|a l | 2 )ia l | 2 . (3) 



We note that the Eq. (2) has more nonlinear terms com- 
pared to the Eq. (1) and, in addition, the Eq. (2) is more 
relevant in condensed matter physics. To the best of our 
knowledge, this equation has not been used to study the 
formation of SL states even though an equation similar 
to Eq. (2) has been used to study the formation SL states 
in one-dimensional system due to a single impurity and 
a symmetric dimer impurity only [ fl8| . Therefore, it is 
essential to consider the right equation (i.e., Eq. (2)) 
for the elaborate study of the formation of SL states in 
one-dimensional systems. Thus we study it here. 

This paper is organized as follows. In Sec. II, we dis- 
cuss the effect of a single impurity (i.e., Xn = X^n.o) on 
the formation of SL states. In Sec. Ill we consider the 
case of a dimer impurity, Xn = {xi S n,o + X2$n,x)- In 
Sec. IV the periodically modulated nonlinear chain as 
well as a perfectly nonlinear chain will be considered. In 
Sec. V, we consider the presence of random nonlinearities 
through out the lattice. The stability of the states are 
discussed in Sec. VI. Finally we summarize our findings 
in Sec. VII. 



II. SINGLE NONLINEAR IMPURITY 

Consider the system of a one-dimensional chain with 
a nonlinear impurity at the cener site. The time evolu- 
tion of an exciton in the system is governed by Eq. (2) 
with Xn = X$n,o- The Hamiltonian which produces the 



1 



equation of motion for the exciton in the system is given 
by 



H - ^(C*C„ + i + C n C* +1 ) 



|(|C 1 | 2 + |C_ 1 | 2 + 2|C | 2 )|Co| 2 . 



(4) 



Since \C n | 2 is a constant of motion, we suitably renor- 
malize so that X)nl^™| 2 = 1- Therefore, \C n \ 2 can be 
considered as the probability for the exciton to be at site 
n. Since we are interested in the stationary localized 
states, we consider the ansatz 



C n = (fr n exp(-iEt); <fr n = (fr V M 



(5) 



where 77 lies between and 1 and can be asymptotically 
defined as r\ = \ E \-^ E ' 2 - ±^ being the energy of the lo- 
calized state which appears outside the host band. Since 
in a one-dimensional system, states appearing outside the 
band at the expanse of the states in the continuum are 
exponentially localized, the ansatz in Eq. (5) is justified 
and may also be readily derived from the Green function 
analysis |t(J[Ilj. Substituting the ansatz to the normal- 
ization condition we get 



1 -rf 
1 



7] 



(6) 



Direct substitution of (fro, <fr n and hence C n in terms of 77 
in Eq. (4) yields an effective Hamiltonian, 



H, 



elf 



4ry x(l-r/ 2 ) 



2 \2 



I + 77 2 (l + V 2 ) 



(7) 



The fixed point solutions of the reduced dynamical sys- 
tem described by H c g will give the values of 77 which 
correspond to the localized state solutions |J. It should 
be noted that since x is constant, the effective Hamil- 
tonian is a function of only one dynamical variable 77. 
Thus fixed point solutions are readily obtained from the 
condition dH c s/drj = 0, i.e., 



— =7?(3 + 7 ? 2 ) = /(7 ? ). 

Ixl 



(8) 



Thus, for a fixed value of x> the number of solutions 
of Eq. (8) lying between and 1 will correspond to the 
number possible SL states. We note that f(rj) is a mono- 
tonically increasing function of 77 and will be maximum 
at 77 = 1. Thus, we see that there is no solution for 
x| < 0.5 and only one solution for \\\ > 0.5. Therefore, 
we conclude that there is a critical value of |x| = 0.5 
separating the no-state region from the one-state region. 
On the other hand, for the same system described by 
Eq. (1), the critical value of |x| separating the no-state 
region from the one-state region is 2. Thus the nearest 
neighbor coupling of the lattice oscillators are responsi- 
ble to reduce the critical value of nonlinear strength by 



a factor of four. Therefore, we expect that next nearest 
neighbor coupling of the lattice oscillators will further 
reduce the critical value of the nonlinear strength and 
thus the critical value may approach to zero when the 
oscillators are globally coupled. 



III. DIMER NONLINEAR IMPURITY 

We consider a one-dimensional lattice with two non- 
linear impurities placed at site and 1 denoted by Xn = 
(Xi^n,o + X2#n,i)- If Xi — Xii we ca ll it symmetric-dimer 
impurity and otherwise asymmetric-dimer impurity. We 
will here consider the system with asymmetric-dimer im- 
purity, since the results for the symmetric-dimer impurity 
will be obtained as a special case. The one-dimensional 
chain with the asymmetric dimer impurity can be de- 
scribed by Eq. (3) with Xn = (xi#n,o + X2<5„,i). For sta- 
tionarity condition we assume that C n — <fr n exp(— iEt). 
Furthermore, for localized states we assume the following 
forms for <fr„ . 



K = [sgn^rT-Vi; 
0_ w = [sgn(£)»j] W to 



n > 1, 
n < 0, 



(9) 



with 77 as defined earlier. The ansatz given by Eq. (9) is 
exact and can be derived from the Green function analy- 
sis. It is also justified from the fact that the states which 
appear outside the host band are exponentially localized 
]10{ . Three different cases are possible: (i) <fr\ = (fro (sym- 
metric), (ii) (fr\ — —<fro (antisymmetric) and (iii) (fri 00 
(asymmetric). To consider all three cases, we introduce a 
variable (3 = 4^. Without any loss of generality, (fro and 
(fri can be assumed to be real. Therefore, the value of (3 is 
confined between 1 and — 1 if \(fro\ < \<f>\\- Else we inverse 
the definition of /3. Thus, (3—1 corresponds to the sym- 
metric states, [3 = — 1 corresponds to the antisymmetric 
states, and — 1 < f3 < 1 corresponds to the asymmetric 
states. Substituting the ansatz in Eq. (9) and the defini- 
tion of f3 in the normalization condition, ~^2 n \C n \ 2 = 1, 
we obtain 



1 -T7 2 

TT^ 2 



(10) 



Using the ansatz and the normalization condition, we 
obtain the effective Hamiltonian as a function of (3 and 
77, given as 



H cS = 2(3 



I-T7 2 

IT^ 2 



2sgn(£')?7 



xi Awr , (i+?? 2 )(i-?7 2 y 



1 + /3 2 



a+/3 2 r 



X2 f ((3 2 +PW + 2(3^(1 - V 2 Y 



(l + f3 2 Y 



(11) 



If j3 = ±1, we obtain 



2 



H± = ±(1 - v 2 ) + 2sgn(E) V 



2 I 2 4 



X2 ( (3 + rP) (l-y 2 ) 
2 I 4 



(12) 



where '+' corresponds to the symmetric states and '— ' 
corresponds to the antisymmetric states. The number 
of fixed point solutions of the reduced dynamical system 
described by H e s gives the possible number of SL states. 
The fixed point solutions satisfy the equation, 



v(i-v 2 Y 



4 (sgn(£) =F ??) 
+ Xi ^(l-rj 2 ) 

+ X2 (v(l-r) 2 ) (3 + /7 2 ) 



+ r?(l-r ? 4 ) 



v(i-v 2 Y 



0. (13) 



one symmetric state and one antisymmetric state on the 
the dotted line. 

Now, let us find the asymmetric states where (5 ^ 1. 
The effective Hamiltonian is a function of two dynamical 
variables, (5 and ry. Therefore the fixed point solutions 
will obey the equations given by 



dH eS 



n , dH cS 
and — — — = 0. 



(14) 



drj d[3 
After a little algebra we obtain the desired equations, 

2 (l-/? 4 )+ Xl /3(l- ?7 2 )(3 + /3 2 + 2r ) 2 ) 

- X2P (1 - V 2 ) (1 + 2/?V + 3/3 2 ) = 0, (15) 
and 

4 fa (1 + /3 2 ) - 2 (1 + /3 2 ) 2 - xiV (1 - (3 + 2/3 2 + 3r? 2 ) 

(16) 



X2V 1 



(2/3 2 + 3/?V + 3/3 4 ) = 0. 
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FIG. 1. Phase diagram of symmetric and antisymmetric 
SL states for the one- dimensional chain with a dimer impu- 
rity. Region I has no SL state, region II has one symmetric 
SL state and region III has one symmetric SL state and two 
antisymmetric SL states. 



After solving this equation for each point in the 
( I Xi I j 1X2 1 ) plane, we obtain the number of possible 
symmetric as well as antisymmetric SL states in the 
(IXi|)IX2|) plane. The phase diagram is shown in 
Fig. 1. The solid line is the critical line separating 
the no-symmetric-state region from the one symmetric- 
state-region and, the dotted line is the critical line 
separating no-antisymmetric-state region from the two- 
antisymmetric-states region. Thus, there is no state in 
region I, one symmetric state in region II, and one sym- 
metric and two antisymmetric states in region III. Fur- 
thermore, there is one symmetric state on the solid and, 




FIG. 2. Phase diagram of asymmetric SL states for the 
one- dimensional chain with a asymmetric-dimer impurity. 
Solid line is the critical line separating no-asymmetric-state 
region (marked as I) from the one-asymmetric-state region 
(marked as II). 



Plotted in Fig. 2 is the phase diagram obtained from the 
solutions of Eqs. (15) and (16). There is only one critical 
line (for |xi| 7^ IX2I) separating the no-asymmctric-statc 
region (marked as I) from the one-asymmetric-state re- 
gion (marked as II). The critical line in Fig. 2, in fact, 
turns out to be very closed to the solid critical line in 
Fig. 1. However, |xi| = | X2 1 is a special far as 

asymmetric SL state is concerned. For |xi| = IX2I, the 
critical value of |x| separating the no-state from the one- 
state region is 2.015. The full phase diagram of SL states 
for the dimer impurity may be obtained by superpos- 
ing Fig. 1 and Fig. 2 (not shown). It is clear from the 
Fig. 1 and Fig. 2 that there are three distinct regions 
having one, two, and four states in the (|xi|, IX2I) plane. 
Furthermore we see that there is a continuity as we go 



3 



from one region to another. For the symmetric dimer 
(lxi| = |X2 1 = |xl)i we obtain three critical points, i.e., 



lY (1) | 
\Xcr | 



Iy (2) I 

\Xcr | 



2.015 and \\ 



(3), 



0.25, 

no SL state for \x\ < |Xer |i one state for 



2.732. There is 



,( 2 ) I 



(2) 



,(3) I 



|y (1) 

\Xcr 



< \x\ < 

three states 



xir } |- The 



|Xcr I, two states tor \xhr \ S |X| < IX 
for |x| = \xcr I an d four states for |x| > 
corresponding critical values obtained by Eq. (1) are 1, 
8/3 and 8 respectively. Thus, in case of symmetric dimer 
also it is observed that the nearest neighbor coupling of 
the lattice oscillators is responsible to reduce the critical 
values, while the maximum number of SL states remain 
unchanged. However, it should be noted that the asym- 
metric nature of the dimer does play a crucial role to 
reduce the critical value of nonlinear strength abruptly. 



IV. FULLY NONLINEAR CHAIN 

We now consider a periodically modulated nonlinear 
chain, i.e., the alternate sites are of the same nonlinear 
strength but any two consecutive sites are of the different 
strength. The distribution of the nonlinear strength, \n 
is given by 



X2n+1 — Xl; X2r, 



X 2 ; n= -oo,...- 1,0,1, ...oo. 



(17) 



When xi — X2, the system reduces to a perfectly nonlin- 
ear chain. In both cases the translational symmetry of 
the system remains preserved. The Hamiltonian for this 
system is given by Eq. (3) with Xn as defined in Eq. (17). 
Using the stationarity condition, C n — <j) n exp(-iEt), we 
can obtain the Hamiltonian in terms of <p n . In this case 
it is not possible to find the exact ansatz for the localized 
states, but there are a few rational choices. For exam- 
ple, the on-site peaked as well as the inter-site peaked 
and dipped solutions are possible. We will consider these 
cases subsequently. 



A. On-site peaked solution 

Since the lattice is periodically modulated with a pe- 
riod of two, the localized solution may have peak either 
at a site having nonlinear strength xi or at a site having 
strength xi- We will carry out the calculation assum- 
ing that the solution is peaked at the center site (site 

0) . The result for the solution peaked at first site (site 

1) may easily be obtained by interchanging the two non- 
linear strengths in the resulting expression for the case 
of the center site peaked solution. Therefore, using the 
ansatz 4> n — 0o^' n ' and the normalization condition we 
get the effective Hamiltonian, 



±^_X2 / V (l-r? 2 ) 2 (1-7? 2 ) (! + ,/») • 
+ V 2 2 I (1 + tj2) 3 (1 + r? 2 ) 3 (1 + 77 4 ) 



V (l- V 2 



4/7 4 (1 - T] 2 



(i + v 2 y (i + ?7 4 ) 



(18) 



After analyzing the fixed point equation, dH c ff/di] = 0, 
we find that there is only one SL state. Interchanging 
Xi with X2i we obtain another SL state. Therefore for 
fixed value of xi and %2, we will have two SL states, 
one peaked at a site of nonlinear strength xi and an- 
other peaked at a site of nonlinear strength X2 ■ Thus we 
conclude that if the nonlinear chain is modulated with 
periodicity n, there will be n number of on-site peaked 
localized solutions. We note that the extra SL state due 
to extra periodicity in the periodically modulated chain 
was overlooked when the same system was studied by the 
use of Eq. (1) 14|. These two states get closer and be- 
come the same as xi ~~ > X2- After substituting Xi~X2^X 
in the fixed point equation we obtain 



1 

Jx\ 



n 



\-J] 2 



(19) 



Equation (19) clearly indicates that there is one on-site 
site peaked solution (SL state) for a perfectly nonlinear 
chain, irrespective of the strength of the nonlinearity. 



B. Inter-site peaked and dipped solution 

For the inter-site peaked and dipped solutions we use 
the dimer ansatz in Eq. (9). After some algebra we ob- 
tain the effective Hamiltonian for the reduced dynamical 
system as 



H cS = 2sgn(S)r ? 



2/3(1 



/3 2 (1-^ 



(1 + /3 2 ) "V(l + /? 2 ) 2 

^(l-Tffll+Z? 4 ) 77 2 (1 - 77 2 ) 2 (1 

(l + (3 2 ) 2 (l-rf) + (i + 02)2(1 - rfi) 
_Xi / 2(l-ry 2 ) 2 (l + /3V) 
2 v (1 + /3 2 ) 2 (1 - r/8) 
2(l-77 2 ) 2 (/? 4 + T7 4 ) 
(l + /? 2 ) 2 (l-7?8) 



X2 

2 



where x—Xi +X2- We first consider the case, (3 
Substituting j3 = ±1 to the Eq. (20) we obtain 

4 77(2 + J] 2 - r) A ) 



|X±| (1T# + jj 2 ) 



(20) 



±1. 



(21) 



from the fixed point equation. Here the '+' corresponds 
to the symmetric SL states and the ' — 'to the antisym- 
metric SL states. From Eq. (21) it is obvious that there 
will always be one symmetric solution for any point in the 
(lxi|i 1X2 1) plane. But for antisymmetric solutions there 
are two critical lines. One being \x2 1 =6.427 — \xi\ ( dot- 
ted line in Fig. 3) and another \x%\ = 8 — \xi\ (solid line 
in Fig. 3). Three distinct regions (in Fig. 3) separated by 
the critical lines are marked by I, II and III. Region I has 
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no SL state, region II has two antisymmetric states and 
region III has only one antisymmetric state. Furthermore 
we note that, for xi = X2 = X; there is no antisymmetric 
state for \\\ < 3.2135, two states for 3.2135 < \x\ < 4, 
and one solution for |%| > 4. These critical values are 
again lower than those obtained from Eq. (1). 



consisting of random nonlinear site energies |19[. The 
evolution of exciton is defined by Eq. (2) with Xn being 
chosen randomly with equal probability for all n. The 
mean square displacement is defined as 

00 

m(t)= J2 * 2 !^WI 2 - ( 22 ) 





lx,l 



10 
t 



FIG. 3. Phase diagram of SL states for a periodically mod- 
ulated nonlinear chain. The regions I, II and III has no, two FIG. 4. The < m(t) > /t 2 is plotted as a function of time 
and one antisymmetric SL state respectively. t. 



Now, we look for the asymmetric solution corresponding 
to (3 7^ 1. The fixed-point equations are lengthy and, 
therefore, are not presented here. The fixed-point equa- 
tions are solved numerically and the possible number of 
asymmetric SL states in the ]X2 1) plane are found. 

It turns out that there is always one asymmetric SL state 
in the (|xi|, | X2 1 ) plane except along the line |xi| = \xz\ 
on which there exists a critical point, |xi| = | X2 1 = 1-75, 
below which there is no asymmetric SL state and above 
which one asymmetric SL state. Thus, it is found that 
the modulation of nonlinearity reduces the critical value 
to obtain a localized state. 

We emphasize that the SL states obtained in the fully 
nonlinear systems arc based on the rational ansatz. How- 
ever, many other SL states may possibly exist in the sys- 
tem. Thus, we are not in a situation to conclude about 
the maximum number of SL states for the fully nonlinear 
chain. 



V. RANDOM NONLINEARITY 

Here we consider a fully nonlinear chain on which all 
sites have random nonlinear strengths. It is known that, 
due to static random site energies, all the states in the 
system are localized and the system behaves as an in- 
sulator. The question is whether or not all the states 
are localized due to random nonlinearities in the system. 
Since the appropriate ansatz is not known, we calculate 
the mean square displacement of a particle in the system 



Since the set of coupled equations [Eq. (2)] can not be 
solved analytically, we solve them numerically using the 
fourth order Runge-Kutta method. Self-expanding lat- 
tice is used to avoid the finite size effect. The particle 
is initially placed at site i.e., the center site. The time 
interval for the calculation is taken to be 0.001. The total 
probability of the particle in the system should be con- 
served and it has been confirmed at every time step of 
our simulation. The mean square displacement is aver- 
aged over hundred random configurations and denoted as 
< m(t) >. Shown in Fig. 4 is < m(t) > /t 2 as a function 
of t. We see that < m(t) > /t 2 initially decreases and 
becomes constant, much smaller than 2. Thus, we con- 
clude that there is a partial ballistic flow of the particle 
and hence all the states are not localized unlike the case 
of a linear system with random site energies. 

VI. STABILITY 

The stability of the SL states can be understood from a 
flow diagram. For this purpose, we consider antisymmet- 
ric solutions due to the presence of a symmetric dimer in 
one dimensional chain. The reduced dynamical Hamilto- 
nian for this case is given by 

H cS = -(1 - v 2 ) + 2r, + M ((3 + v 2 ){ i _ ,2)2). (23 ) 

The flow diagram of the dynamical system described by 
the iJ c ff is constructed in the following manner. We treat 
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dH c g/drj — f(rj) as the velocity and rj as the coordinate 
of the dynamical system. We fix \ to be 3 for which two 
SL states appear. f(rj) is plotted in Fig. 5 as a function 
of -q. A and B are the fixed points corresponding to the 
SL states. If f(rj) > 0, the flow of the dynamical variable 
is in the right direction, or else, in the left direction, as 
indicated by arrows in the figure. It is clear from the 
flow diagram that A is a stable fixed point, whereas B is 
unstable. Therefore, of the two antisymmetric SL states 
one is stable and another one is unstable. Similarly, we 
analyze the stability for all the states. It is found that if 
N states exist for a particular case, N/2 states arc stable 
provided that N is even and, otherwise (N+l)/2 states 
are stable. 



0.5 








n 



FIG. 5. The flow diagram where f{rf) is plotted as a func- 
tion of rj. A and B are the fixed points. Point A is a stable 
while B is unstable. 



VII. CONCLUSION 

Modified discrete nonlinear Schrodingcr equation has 
been used to study the formation of stationary localized 
states in a one-dimensional system with a single nonlin- 
ear impurity and a dimer (symmetric as well as asym- 
metric) nonlinear impurity. It is found that the critical 
values of nonlinear strength decrease due to the nearest 
neighbor coupling of the lattice oscillators. Furthermore, 
a periodically modulated nonlinear chain is studied and 
new SL states are found due to the modulation of non- 
linearity. It is important to note that there will be more 
number of SL states if the periodicity of the modulation 
is larger. Another interesting point to note is that the 
critical value of nonlinearity to obtain SL state decreases 
drastically due to the modulation of the nonlinearity in 
a fully nonlinear chain and the asymmetric nature of the 
dimer impurity. The phase diagrams of SL states are pre- 
sented for all the cases. The mean square displacement 
of a particle in one-dimensional chain with random non- 
linear site energies is calculated to investigate whether 



all the states are localized as it happens in the case of a 
linear system. However, it is found that a fraction of the 
particle flows ballistically indicating that all the states 
are not localized even in the presence of random nonlin- 
earities through out the system. It is further found that 
SL states may appear even in a perfectly nonlinear sys- 
tem. Thus, one may call these SL states as self-localized 
states. The stability of the SL states is also discussed. To 
see the effect of these nonlinearities on the formation of 
SL states in lattices with higher connectivity and higher 
dimension, the Cayley tree may be a good candidate as 
a linear host lattice. 
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